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Abstract 

We investigate the influence of oxygen annealing on the room temperature optical 
and scintillation properties of CaW04 single crystals that are being produced for 
direct Dark Matter search experiments. The applied annealing procedure reduces the 
absorption coefficient at the peak position of the scintillation spectrum (~ 430 nm) 
by a factor of ~ 6 and leads to an even larger reduction of the scattering coefficient. 
Furthermore, the annealing has no significant influence on the intrinsic light yield. 
An additional absorption occurring at ~ 400 nm suggests the formation of O hole 
centers. Light-yield measurements at room temperature where one crystal surface 
was mechanically roughened showed an increase of the measured light yield by ~ 
40% and an improvement of the energy resolution at 59.5 keV by ~ 12% for the 
annealed crystal. We ascribe this result to the reduction of the absorption coefficient 
while the surface roughening is needed to compensate for the also observed reduction 
of the scattering coefficient after annealing. 

1 Introduction 

Calcium tungstate (CaW04) was one of the first known scintillators. Its radiolumines- 
cence was discovered by Thomas Edison in 1896 pp. However, it is not widely used in 
particle physics mainly because of its long luminescence decay times in the order of us 
[2]. Nevertheless, CaW04 is an interesting target material for direct Dark Matter search 
experiments because of its relatively low intrinsic radioactivity and high light yield. In 
addition, the heavy atomic mass of tungsten (A = 183) makes it a favorable target 
for coherent (spin-independent) WIMP-nucleus scattering^] because of the expected A 2 - 
dependence of the cross section [3J. 



1 WIMP: Weakly Interacting Massive Particle 
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The CRESST-II (Cryogenic Rare Event Search with Superconducting Thermometers) 
experiment attempts to detect WIMPs using CaWC-4 crystals operated as low-temperature 
detectors [3J. With these detectors an active background discrimination is achieved by 
the simultaneous detection of phonons and scintillation light produced by particle inter- 
actions in the crystal. 

The most crucial part in the background discrimination ability of a CRESST-II detector 
module is the amount of detected scintillation light. Therefore, the crystals should have 
a high intrinsic light yield. The intrinsic light yield is defined as the number of scin- 
tillation photons that are produced per MeV of deposited energy. Since the cylindrical 
crystals have a large size of ~ 40 mm diameter and ~ 40 mm height also scattering and 
absorption of scintillation light in the crystal has to be taken into account. This is of 
particular importance as in a CRESST-II detector module, there is no optical coupling 
between the scintillating crystal and the light detector. Due to the high symmetry of 
the cylindrical shape and the high refractive index of CaW04 (n~ 1.95) ~ 50% of all 
photons will always be internally reflected [5] and are thus trapped inside the crystal. 
This means that the measured light yield which is defined as the number of detected 
photons for a given energy is usually much lower than the intrinsic light yield. Simula- 
tions have shown that scattering reduces the fraction of trapped light [6j since a random 
change of the direction prevents a photon from being always internally reflected. In this 
prospect, the crystals should exhibit a large scattering coefficient and a small absorption 
coefficient. 

The EURECA (European Rare Event Calorimeter Array) project is a future ton-scale 
experiment which combines efforts of cryogenic Dark Matter search in Europe (CRESST, 
EDELWEISS^, ROSEBUD^) [7j. It is planned that a large part of the target mass will 
consist of CRESST-like detector modules. 

For future runs of CRESST and for the EURECA experiment it is crucially important 
to ensure the availability of CaW04 crystals, which fulfill the requirements regarding 
their optical and scintillation properties as well as their radiopurity. It is therefore desir- 
able to have a direct influence on the selection of the raw materials (CaC03, WO3), the 
crystal growth and the after-growth treatment. For this reason, a Czochralski furnace 
dedicated to the growth of CaW04 single crystals has recently been installed at the 
Crystal Laboratory of the Technische Universitat Miinchen (TUM) |8j. 
It is known that the thermal annealing under oxygen atmosphere can ameliorate the 
mechanical, optical and luminescent properties of CaW04 crystals [TDJ [TT]. Oxy- 
gen annealing is also used to improve the properties of other tungstate scintillators like 
ZnWC>4 and PbW04 |T21 E3] . Especially PbW04 was extensively studied since it is used 
for the electromagnetic calorimeter of the CMS (Compact Muon Solenoid) experiment 
[14j . High-temperature annealing results in relief of internal stresses. In addition, oxy- 
gen annealing can reduce the oxygen deficiency which is present after growth due to the 
high growth temperature (~ 1600 °C) and the reduced oxygen partial pressure in the 
growth atmosphere (see Section I27T]) . 

2 Experience pour DEtecter Les Wimps En Site Souterrain 
3 Rare Objects SEarch with Bolometers UndergrounD 
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In this work we use a combination of experiments and Monte-Carlo simulations to de- 
termine the absorption and scattering coefficients as well as the intrinsic light yield of 
a CaWC>4 crystal at room temperature before and after oxygen annealing. In addi- 
tion, light-yield measurements at room temperature were conducted to determine the 
influence of the annealing process on the performance of the crystal as a scintillator. 



2 Experimental 

2.1 Sample Preparation 

For the measurements a cube-shaped crystal of 18 x 18 x 18 mm 3 size was used which was 
produced from an ingot grown by the Czochralski method at the Crystal Laboratory of 
the TUM. The ingot was grown in a rhodium crucible under constant flow of a mixture 
of 99 % argon and 1 % oxygen. The details of the growth process are described elsewhere 
[HI [15] . All surfaces of the crystal were optically polished. The annealing was carried out 
at a temperature of 1450 °C under constant flow of pure oxygen and lasted 48 h. The 
crystal was heated up and cooled down at a rate of 100°C/h. 



2.2 Transmission Measurement 



For the transmission measurements a Perkin Elmer LAMBDA 850 UV/VIS spectropho- 
tometer was used. The transmittance was measured for wavelengths from 250-800 nm 
in steps of 1 nm. The transmittance T is defined as 



T = h/h 



(1) 



where I\ and Iq are the measured intensities with and without the sample in the beam, 
respectively. The attenuation coefficient a a tt was obtained using the following equation 



(1 - R) 2 ■ exp(-a att d) 
1 - R 2 ■ exp(-2a aU d) 



Oiatt 



In 



-(1-R) 2 + ^(1 - fi) 4 + 4T^R2 
2TR 2 
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(2) 



(3) 



Here d is the thickness of the crystal and R its reflectivity. The denominator of Equation 
([2]) accounts for multiple reflections. 

CaW04 (space group: I^i/a) is weakly birefringent (5 ss 0.017 |16j), however, in the 
measurements the crystal was aligned with the beam of light parallel to the optic axis 
(c-axis) of the crystal so that no birefringence occurred. Therefore, the reflectivity R 
can be calculated as 

(n - l) 2 



R 



K + 1) 2 



(4) 
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with n being the ordinary refractive index which was calculated by the following dis- 
persion formula [17j : 

2 _ 2.5493 • (A//im) 2 0.92 • (A//xm) 2 

U °~ ~ {X/fim) 2 - 0.1347 2 + (X/nm) 2 - 10.815 2 ^ ' 

where A is the photon wavelength. 



2.3 MCRIM Technique 

The Monte-Carlo Refractive Index Matching (MCRIM) technique is a combination of 
light-yield measurements and Monte-Carlo (MC) simulations to determine unknown op- 
tical and scintillation properties of heavy inorganic scintillators at room temperature 
|18j . For a crystal with known attenuation coefficient a a a it can be used to determine 
the intrinsic light yield Lq as well as the ratio B of the scattering coefficient a sca t and 
the absorption coefficient a abs '- 

B = a scat /a abs (6) 

OLatt = a abs + a scat (7) 

Using Equations ([6]) and ([7]) one can then calculate the values for the scattering and 
absorption coefficients. 

The MCRIM technique is based on the principle that the measured light yield L m of 
a scintillating crystal is defined as the product of the intrinsic light yield Lq and the 
light collection efficiency ij of the setup in which the crystal is placed. In general, r\ 
is dependent on a number of parameters such as the geometry of the crystal and the 
setup, as well as the ratio B and other optical properties of the materials used in the 
experiment: 

L m = r]L (8) 

As the measured light yield is dependent on both B and Lq, these quantities cannot be 
determined in a single experiment. However, by taking the ratio of the measured 
light yields L m ^\ and L m ^ from two distinctive measurements with different collection 
efficiencies rji and 772 the dependence on Lq can be eliminated: 

r 1/2 = L m ,i/L m ,2 = m/m (9) 

The experimental setup of the MCRIM technique (see Fig. [1]) comprises the crystal 
the properties of which are to be determined and a photomultiplier tube (PMT) placed 
together in a light-tight box. A small gap between the crystal and the window of the 
PMT offers the possibility of introducing materials of different refractive indices n. The 
ratio T\j2 is obtained by measuring the light yield with two different materials in the gap 
when the crystal is irradiated with a radioactive source. MC simulations of the setup 
that iterate the value of B can then be used to match the measured ratio r\j2 if all 
other optical properties are known. In this prospect, the experimental setup is designed 
to minimize the presence of optical components such as unpolished surfaces and non- 
ideal reflectors, which are difficult to implement accurately in a simulation. Having 



4 



determined B through MC matching, the collection efficiency r\ can be calculated by the 
MC. Using Equation ([8]) it is then possible to deduce the intrinsic light yield Lq. A more 
detailed description of the MCRIM technique is given in Ref. [18J. The MC simulations 
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Fig. 1: Experimental setup of the MCRIM technique. The crystal is placed onto a photomultiplier 
tube (PMT) and irradiated with a gamma source. A small gap between crystal and PMT offers 
the possibility to introduce materials with different refractive indices. 

of the setup were performed using the GEANT4 (version 4.9.2p01) software [19]. For the 
measurements an ETL 9305KB photomultiplier was used which was operated in photon 
counting mode. The PMT pulses were recorded by an Acqiris DC282 digitizer with a 
sampling rate of 1 GHz and analyzed off-line. The crystal was irradiated by a 241 Am 
gamma source (59.5 keV). At first a measurement was performed just with air (n a j r = 1) 
between the PMT and the crystal. For the second measurement the gap was filled with 
an optical geQ (n ge i = 1.46). The ratio r air / ge i is defined as the measured light yield 
L m ,air with air divided by the light yield L m ^ ge i measured with the gap filled with the 
gel: 



2.4 Light- Yield Measurements 

To investigate the influence of the annealing procedure on the performance of the crystal 
as a scintillator additional light-yield measurements were performed at room tempera- 
ture. In these measurements the crystal was mounted in an aluminum housing (see Fig. 
[2]) which was then placed onto the PMT of the setup described in Section T2.3I The hous- 
ing is covered with a highly reflective foiH on the inside and the crystal is supported by 
small holders made of polyethylene. In this configuration there is a small gap between 
the crystal and the side walls of the reflective housing as well as between the crystal 
and the PMT window (see Fig. [2|). In this way the geometry is similar to that of a 

4 Thoriabs G608N refractive index matching gel 
5 3M Radiant Mirror Film VM2000 
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CRESST-II detector module 



It is known that surface roughening can increase the 
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Fig. 2: Experimental setup of the additional light-yield measurements. The crystal was mounted 
in a reflective housing and placed onto the PMT window in a geometry which is similar to that 
of a CRESST-II detector module. 



light collection especially when there is no optical coupling between crystal and light 
detector [20]. This increased light collection can be ascribed to the fact that each pho- 
ton will hit the roughened surface at a random incident angle which prevents it from 
being always internally reflected. Therefore surface roughening decreases the fraction 
of light trapped in the crystal similar to scattering. Furthermore, a rough surface can 
show strong distributions of transmitted light for angles of incidence larger than the 
critical angle where no transmitted light would be observed at a specular surface |21| . 
In CRESST-II, a significant improvement of the energy resolution was obtained by me- 
chanical roughening of the crystal surface that is facing the light detector using silicon 
carbide powder with a grain size of ~ 9 pm [22]. Therefore, in the measurements that 
are described in this section the crystal surface that is facing the PMT was roughened 
by applying the identical procedure]! The crystal was again irradiated with a 241 Am 
gamma source (59.5 keV). 



3 Results and Discussion 

3.1 Transmission Measurement 

Fig. [3] shows the measured transmittance T (see Equation [1]) of the crystal before and 
after annealing in dependence of the photon wavelength A. It can be seen that the 
transmittance is considerably increased after the annealing process. Fig. [5] shows the 
attenuation coefficient a a tt calculated from the measured transmittance (see Equation 
ED versus the photon wavelength. Tab. [T] summarizes the values of the attenuation 

6 The roughening procedure was also reapplied after annealing. 
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coefficient at a photon wavelength of 430 nm which corresponds to the approximate 
peak position of the CaWC>4 scintillation spectrum at room temperature [23J. 
The error for the attenuation coefficient is given at 95 % CL and was calculated from the 
data of four independent transmission measurements. However, it has to be pointed out 
that, in addition, there are some systematic uncertainties in the value of a a tt because 
the refractive index of the crystal was not directly measured. Furthermore, surface 
irregularities can cause errors in the measured transmittance in the order of ~ 20% 
[6]. Tab. [1] shows that the attenuation coefficient at 430 nm is decreased by a factor 
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Fig. 3: Transmittance T of the CaWC>4 crystal before and after annealing versus photon wave- 
length A. The pronounced absorption around 400 nm in the annealed crystal is probably caused 
by O - hole centers. 





a a tt [cm x ] at 430 nm 


before annealing 
after annealing 


0.335 ± 0.026 
0.039 ± 0.010 



Tab. 1: Value of the attenuation coefficient a a tt at 430 nm of the CaWC>4 crystal before and after 
annealing. The wavelength corresponds to the approximate peak position of the scintillation 
spectrum at room temperature. The errors are given at 95 % CL and were determined from the 
data of four independent transmission measurements. 



of ~ 8 after the annealing procedure. This can be attributed to the removal of crystal 
defects which otherwise cause scattering or absorption of the photons. Scattering centers 
in CaWC>4 crystals include solid and gaseous inclusions [251 E5]- Solid inclusions can 
be caused by the incorporation of rhodium from the crucible [26} 127] while gaseous 
inclusions probably result from oxygen separation from WO4 - complexes [28]. It is, 
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Fig. 4: Attenuation coefficient a a tt of the CaWC>4 crystal before and after annealing versus 
photon wavelength A. 



however, unclear if inclusions that are present in the crystal can be removed by annealing. 
Dislocations that are present in the crystal can act as dielectric inhomogeneities and 
hence also scatter light [29]. The formation of dislocations can result from temperature 
gradients after growth in the cooling crystal [25]. During annealing dislocations become 
mobile and can be removed from the crystal or form low-angle boundaries [30} [25] . This 
may be a reason for the observed reduction of the attenuation coefficient. In addition, 
the main intrinsic defect type in CaWC>4 crystals are oxygen vacancies [9] and their 
formation is facilitated by the reduced oxygen partial pressure in the growth atmosphere 
(see Section I2.ip . These oxygen vacancies can be filled during annealing under oxygen 
atmosphere which can also reduce optical scattering and absorption. 
In Fig. [3] one can observe a pronounced absorption around 400 nm in the annealed 
crystal. According to Ref. [ID] oxygen can also be incorporated on interstitial sites 
during annealing which leads to absorption bands at ~ 400 nm and ~ 310 nm. However, 
it is also known that for PbWC>4 an absorption band at ~ 420 nm can be induced after 
oxygen annealing [3HI32j. This absorption band is commonly ascribed to O - hole centers 
|31j . The filling of oxygen vacancies during annealing might lead to the creation of O 
centers to maintain local charge balance. It is known that for CaWC>4 the localization 
of holes at one oxygen of the WO4 tetrahedron is the dominant trapping mechanism 
[33]. In this way the O - hole center may be a more natural explanation for the observed 
absorption band. 

The vanishing transmittance below ~ 300 nm in both curves of Fig. [3] is due to intrinsic 
absorption corresponding to the excitation of electrons from the valence band to the 
conduction band of CaWC>4 |23j. 
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3.2 MCRIM Technique 

Tab. [2] shows the results of the light-yield measurements from the MCRIM technique 
before and after annealing of the crystal. The results for the intrinsic light yield Lo, 
the absorption coefficient a a b s and the scattering coefficient a sca t, that were determined 
from the matching of the values of Tab. [2] to MC simulations, are summarized in Tab. El 
All errors are given at 95 % CL and were determined from the data of four independent 
light-yield measurements and the uncertainty of the attenuation coefficient (see Tab. [T]). 
The values show that the annealing procedure leads to a decrease of the absorption 





L m ,air [p.e.] at 59.5 keV 


L m ,gei [p-e.] at 59.5 keV 


before annealing 
after annealing 


39.9 ±0.7 
39.2 ±0.4 


79.3 ±1.7 
92.0 ± 1.0 



Tab. 2: Results of the light-yield measurements from the MCRIM technique performed with a 
241 Am gamma source (59.5 keV). L mMi i r and L m ge i denote the measured light yields in pho- 
toelectrons (p.e.) with the gap between crystal and PMT (see Fig. [1} filled with air and gel, 
respectively. All errors are given at 95 % CL and were determined from the data of four inde- 
pendent light-yield measurements. 





L [ph/MeV] at 59.5 keV 


otabs [cm" 1 ] 


a scat [cm" 1 ] 


before annealing 
after annealing 


24800 ± 3300 
20600 ± 900 


0.231 ±0.051 
0.036 ±0.010 


0.104 ±0.048 
0.004 ± 0.003 



Tab. 3: Results of the MCRIM technique for the intrinsic light yield Lq at 59.5 keV, the absorption 
coefficient a a b s and the scattering coefficient a sca t . All errors are given at 95 % CL and were 
determined from the data of four independent light-yield measurements and the uncertainty of 
the attenuation coefficient (see Tab. [I]). 

coefficient of the crystal by a factor of ~ 6 and an even larger reduction of the scattering 
coefficient. The measurements also suggest that the attenuation of scintillation light in 
the annealed crystal is dominated by absorption. However, as the errors of the scattering 
and absorption coefficients are rather large no definitive statements can be made. 
As mentioned in Section [1] scattering can increase the light collection due to the reduction 
of light being trapped inside the crystal. This explains the fact that the measured light 
yield L TO>a j r with no optical coupling between crystal and PMT is not increased after 
annealing despite the reduction of the absorption coefficient. In contrast to this, the 
measured light yield L m ^ ge i with optical coupling - where trapping is far less prominent - 
increases, as it is hardly influenced by the reduction of the scattering coefficient. Since in 
a CRESST-II detector module there is also no optical coupling to the light detector (see 
Section I2.4[) it would be desirable to increase the transmittance of the crystal without 
removing the scattering centers. However, the reduction of trapped light can also be 
accomplished by surface roughening (see Sections 12.41 and 13. 3|) . 

The value of the intrinsic light yield Lq shows a small, although not significant, decrease 
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after the annealing procedure. This may result from the formation of defects which act 
as quenching centers, i.e. they are responsible for energy absorption followed by non- 
radiative decay. 

We note that the values of the intrinsic light yield from Tab. [3] are comparable to the 
value of 22700 ± HOOph/MeV which was determined in Ref. [18J for a CaW04 crystal 
produced by the Institute of Materials SRC "Carat" (Ukraine) 

3.3 Light-Yield Measurements 

Fig. [5] shows two 241 Am spectra measured with the CaW04 crystal before and after 
annealing. As mentioned in Section 12.41 the crystal surface that was facing the PMT 
had been mechanically roughened for these measurements as it is done for CRESST-II 
detectors. 

In Tab. 0]we show the mean values and errors (95 % CL) for the measured light yield Lm 
as determined from the peak position, and for the energy resolution as derived from the 
ratio of the peak width (FWHM) and peak position. The results were calculated from 
the data of four independent light-yield measurements. In contrast to the measured 
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Fig. 5: Spectra of 241 Am gamma rays (59.5 keV) measured with the roughened CaWO*4 crystal 
before and after annealing. The measured light yield (peak position) and energy resolution 
(FWHM/peak position) were determined from a Gaussian fit to the photopeak. 

light yield L m>a i r of the polished crystal (see Tab. [2|), the measured light yield L m of 
the roughened crystal (see Tab. 0|) is increased by ~ 40% and the energy resolution 
is improved by ~ 12 % after annealing despite of a large reduction of the scattering 

7 The " Carat" institute is one of the current suppliers for CRESST-II detector crystals. 




10 





L m [p.e.] at 59.5 keV 


resolution [%] at 59.5 keV 


before annealing 
after annealing 


31.0 ±0.6 
43.3 ±0.5 


52.2 ±0.6 
46.0 ± 1.0 



Tab. 4: Measured light yield L m in photoelectrons (p.e.) and energy resolution at 59.5 keV of 
the roughened crystal before and after annealing. The mean values and errors (95 % CL) were 
determined from the data of four independent light-yield measurements. 



coefficient (see Tab. [3])|j As mentioned in Section 12.41 surface roughening can reduce 
the fraction of scintillation light that is trapped in the crystal similar to scattering. We 
thus infer that the effects of the roughened surface are dominant over those caused by 
scattering in the crystal. In this way, the reduction of the scattering coefficient can be 
compensated by the roughened surface and the observed improvement of the measured 
light yield and energy resolution after annealing can be explained by the reduction of 
the absorption coefficient (see Tab. [3]). To support this argument we have used the 
MC simulations mentioned in Section 12.31 to determine the fraction of trapped light in 
the crystal before and after annealing. This was done for the cases that all surfaces of 
the crystal are polished and that one of the surfaces is roughened. In the simulation a 
rough surface is a collection of micro-facets with a Gaussian distribution of facet slopes 
|34j . The standard deviation of this distribution that was implemented in the simulation 
was determined from the measured surface profile of a roughened crystal. The results 





all surfaces polished 


one surface roughened 


before annealing 
after annealing 


62.8 ± 5.2% 
56.8 ±3.2% 


59.8 ± 4.3% 

35.9 ± 4.9% 



Tab. 5: Simulation of the fraction of light that is trapped inside the crystal before and after 
annealing. The simulation was conducted with either all surfaces polished or with one surface 
roughened. The mean values and errors (95% CL) were determined from the uncertainties on 
the absorption and scattering coefficients. 

of the simulations are shown in Tab. It can be seen that the polished crystal shows 
no significant change of the fraction of trapped light after annealing. However, with one 
surface roughened the fraction of trapped light is reduced by ~ 40% after annealing. 
This agrees well with the measured increase of L m by ~ 40% after annealing of the 
roughened crystal (see Tab. Therefore the simulations support the interpretation 
that the effects of surface roughening are dominant and compensate the decrease of the 
scattering coefficient. 

We note that a similar result was found for ZnWO"4 crystals which have shown an 

8 The fact that the measured light yield L m of the roughened crystal before annealing is lower than 
the value L m ,air of the polished crystal (see Tab. [2} results from the different geometries of both setups 
(see Figs. [T]and[2]). In particular, the additional holders of the crystal and the smaller surface area 
where photons can escape the reflective housing and be detected lead to a lower light collection. 
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improvement of the measured light yield in a similar setup by ~ 30% after oxygen 
annealing [13) . 

4 Summary and Conclusion 

We have investigated the room temperature optical and scintillation properties of an 
18 x 18 x 18 mm 3 CaW04 crystal that was produced by the Czochralski method at the 
Crystal Laboratory of the TUM in a setup which is dedicated to the growth of CaW04 
crystals for the direct Dark Matter search experiments CRESST and EURECA. Using 
transmission measurements and the MCRIM technique we determined the absorption 
and scattering coefficients as well as the intrinsic light yield of the crystal. The mea- 
surements were performed before and after the crystal was annealed under a constant 
flow of pure oxygen at 1450 °C for 48 h. 

It was shown that annealing decreases the absorption coefficient of the crystal by a 
factor of ~ 6 and leads to an even larger reduction of the scattering coefficient. This 
result is likely to be explained by the removal of dislocations and oxygen vacancies that 
are present in the crystal directly after growth. Furthermore, the applied annealing 
procedure has no significant influence on the intrinsic light yield. An additional absorp- 
tion occurring around 400 nm in the transmission curve suggests the formation of O 
hole centers during annealing. Further light-yield measurements at room temperature 
were performed where one crystal surface was mechanically roughened as it is done for 
CRESST-II detectors. The results showed an increase of the measured light yield by 
~ 40% and an improvement of the energy resolution at 59.5 keV by ~ 12% for the 
annealed crystal. We ascribe this result to the reduction of the absorption coefficient 
while the surface roughening is needed to compensate for the also observed reduction of 
the scattering coefficient after annealing. 

Our measurements have shown that the oxygen annealing of CaW04 crystals can lead 
to a significant improvement of their performance as a scintillator at room temperature. 
This behaviour has also to be confirmed at low temperatures (~10mK) in the operating 
range of CRESST/EURECA detector modules. 
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